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fi ' The projected SO(5) (pSO(5)) Hamiltonian incorporates the quantum spin and supercon- 

ducting fluctuations of underdoped cuprates in terms of four bosons moving on a coarse 
grained lattice. A simple mean field approximation can explain some key feautures of 
the experimental phase diagram: (i) The Mott transition between antiferromagnet and 
superconductor, (ii) The increase of T c and superfluid stiffness with hole concentration x 
and (iii) The increase of antifcrromagnetic resonance energy as yjx — x c in the supercon- 
f*0 ' ducting phase. We apply this theory to explain the "two gaps" problem in underdoped 

■^4- | cuprate SNS junctions. In particular we explain the sharp subgap Andreev peaks of 

^\ , the differential resistance, as signatures of the antifcrromagnetic resonance (the magnon 

r™^ ■ mass gap). A critical test of this theory is proposed. The tunneling charge, as measured 

(«»»s ■ by shot noise, should change by increments of AQ = 2e at the Andreev peaks, rather 

than by AQ = e as in conventional superconductors. 
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1. Introduction 

The underdoped regime of high T c cuprate superconductors exhibits a pairing energy 
(called pseudogap A p ) which is much larger than the superconducting- transition 
temperature T c . A p is measured by NMFM, tunnelinga, photoemmissionH and other 
probes, as a sharp decrease in the single particle density of states below A p at 
temperatures below 2A p . Assigning A p to the local pairing energy, is consistent 
with the observed anisotropy of |A p (k)| in momentum spaceO as expected for a 
d-wave symmetry of the pair wavefunction. 

In sharp -contrast to conventional BCS superconductors, T c and A p are not 
proportionalQ. While A p decreases with hole doping x, T c grows with x, as does the 
superfluid density p s . The latter relation is not coincidental as argued by Emery and 
KivelsonO. Basically, since p s is small at low doping, superconductivity is destroyed 
by thermal phase fluctuations, and thus T c is determined by the Kosterlitz-Thouless 
transition temperature T c ~ p s , and not by pair breaking effects which are supressed 
at T c « 2A p . 

The behavior of p s ~ x is even more unconventional. It demonstrates that 
the correct number of charge 2e bosons is proportional to the number of doped 
hole pairs away from the half filled Mott phase. The Mott phase exhibits long 
range antiferromagnetic (AFM) correlations. This is a consequence of the strong 
local Hubbard interactions, which cannot be easily handled as a weak coupling 
perturbation of the electron gas. 

The proximity of the AFM phase to the superconductor suggests that the ef- 
fective Hamiltonian for the low energy collective modes would describe 4 bosons: 
three magnons of the Heisenberg antiferromagnet, and the hole pairs. In fact, 
in addition there are gapless fermionic particle-hole excitations near the nodes of 
A P (k),kw(±7r/2,±7r/2). 

In the AFM phase, one magnon branch condenses, yielding a finite staggered 
magnetization, while the other two become the two spin wave modes. At higher 
doping, in the superconducting phase, AFM order disappears which allows these 
modes to become massive with mass A s , and for optimally doped YBCO they are 
observed by inelastic neutron scattering as a resonance at around 40meV near the 
AFM wavevector. r. 

Such as model was introduced!! as the "projected SO (5) (pSO(5)) theory". This 
model is constructed by projecting out doubly occupied configurations from the 
local two-site SO (5) multiplets. The remaining states are local singlets, triplets and 
hole pairs, which define the local bosons. 

In the following, we review the pSO(5) model and its primary consequences 
by mean field theory. The important predictions are p s ~ x, and that A s (x) ~ 
y/x — x c . 

We then use the pSO(5) model to predict current singularities in Superconductor- 
Normal-Superconductor (SNS) junctions!], and propose a future experimental test 
of this theory. The tunneling charges, as measured by shot noise, should change by 
2e at the Andreev peaks, rather than by e as in conventional superconductors. 

2. The Model 

The large onsite Hubbard repulsion between electrons is imposed by an apriori 
projection of doubly occupied states from the Hilbert spaced. 

The undoped vacuum | V) is a half filled Mott insulator in a quantum spin liquid 
state. The pSO(5) vacuum possesses short range antiferromagnetic correlations. A 
translationally invariant realization of \V) on the .microscopic square lattice, is the 
short range resonating valence bonds (RVB) stateE30. Since local singlets without 
holes and double occupation can only be produced with even number of sites, we 
must choose a coarse grained lattice in order to define the vacuum and the Fock 
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space. 

If one defines a superlattice of dimers on the original square lattice, a Fock 
vacuum can be explicitly constructed in terms of electronic operators as dimer 
singlets, where i labels a dimer. 

\V)i = ^= (4, T 4,| - 4 U 4, T ) I0)i (!) 

where 1 0)^ is an empty state of electrons. The hole pairs are simply 

bl\V)i= |0)< (2) 

The triplets are defined as bosonic excitations of the RVB vacuum i.e. 

b l\vh = ^(4, t 4,i + c !u4,t) |o)< 
bUvh = ^(4, T 4,T + 4a4,i) \o)t 
b U v )i = ^(4 )T 4,i-4,i4,t) \°)i ( 3 ) 

Out of the undoped RVB vacuum |V), b h create charge 2e bosons (hole pairs) 

and b' a , a — x,y,z create a triplet of antiferromagnetic, spin one magnons. For 
the square lattice, a similar construction which preserves the lattice rotation group 
and has explicit d-wave|-symmetry of the hole pairs bosons, will be presented in a 
forthcoming publicationu. 

The lattice pSO(5) Hamiltonian is 



ojpSO(5) njcharge , njspin _■_ n/int _j_ n/Coul , *)-/ f ■ 
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w** = wj2--Wht+J2 b * b «i 



(4) 



where : () : denotes normal ordering, and nf = (b ia + b ia )j\[2 is the Neel spin field. 
7i mt describes short range interactions between bosons, and ?{ Coul describes the 
long range Coulomb interactions. ij/ e ™ describes coupling to the nodal (fermionic) 
quasiparticles, which contribute a finite density of single electron states at low 
energies. We shall not compute the fermion contributions, and will discuss them in 
detail elsewhereQ. 

3. Mean Field Theory: Results 

The uniform mean field approximation to Eq. (|4|) is straightforward^. It amounts 
to replacing 6 7 i — > (6 7 ), 7 = h,ot. The order parameters are related to experi- 
mental observables: the Bose condensate of hole pairs is the superconducting order 
parameter 
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x = ( b i) = (J2 d v c h c u) ( 5 ) 

ij 

where dij is the normalized short range pair wavefunction with d wave symmetry. 
Bose condensation of magnons yields the staggered magnetization 

V = %) = (n z ) (6) 

The T = variational energy of H charge + H spm + H mt is 

e mft /m = (£c _ 2 ^_ z_ Jc)x 2 + (£g _ zJs)y 2 + w{x 2 + y2)j (?) 

where z = 4 is the square lattice coordination and A/" is the lattice size. Minimizing 
E(x, y) we find a first order transition between two phases: 

fj, < [i c AFM insulator: i = 0,!/^0 

fi > [i c d-SC: x 7^ 0,y = 



where 



^ ((e ~ e«) ~ (^c - zJ.j) (8) 

At /i < /i c we have an undoped Mott insulator with no hole pair bosons, and 
where the magnons Bose-condense. The condensate supports a finite staggered 
magnetization 

\(n a )\ 2 = (2J s -±e s )/W = m 2 s p, < Mc (9) 

Expanding the mean field Hamiltonian to second order in the Bose operators, one 
obtains two linear spin wave modes at 



qq| 



2V2J s /h. (to) 



c is the semiclassical spinwave velocity which agrees with semiclassical limit of the 
Heisenberg antiferromagnet. 

At fi > [i c the ground state becomes doped with hole pairs which Bose-condense 
into a superconducting phase with an order parameter 



\{bl)\=(»-»c)/W + mi m>Mc (11) 



The mean field phase stiffness is given by p c = J c {b\) 2 , and therefore Eq. ( |ll| ) 
explains why p c increases with chemical potential (,aiid doping) in the underdoped 
superconducting regime, as observed experimentallyEx Quantum phase fluctuations 

which increase with p c significantly reduce the mean field order parameter ( pi] ) 
near the transition. 

Long range interactions in Ji Coul ; frustrate the first order transition and create 
intermediate (possibly incommensurate) phaseso, which we shall not discuss here. 

Analysis of the linear quantum fluctuations about mean field theory in the su- 
perconducting phase, when y — 0, yields three massive magnons (of spin 1). Their 
mean field Greens function approximates the spin structure factor 

S aa >(u,q) ra S0-5 aa t-j (12) 
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Here c is the spin wave velocity, and so is a normalization factor. The poles of Eq. 
( |l2| ) have a mass gap at A s . The mean field magnon gap is found to depend on the 
chemical potential as follows: 



A s 
c(//) 



2a/(/x - fJ.c)(^ - Mc + 2J S ) ex V^ 7 " 
2V2J,Vl + (M-Mc)/(2J s ) 



(13) 



which by Eq. (O) also implies that A 2 , increases, and the magnon dispersion stiffens 
at higher doping. 

Thus the pSO(5) mean field theory can explain, the systematic increase of A s 
with T c which has been obsep/ed by Fong el. alt3. The doping dependent reso- 
nance energy A s (<5) increasesliij between A s (0.5) = 25meV, (with T c = b2°K), and 
A s (l) = 40meV, (at T c = 92°K)B. 



4. SNS Junction: KBT theory 

Peaks in the differetial resistance of Superconducting-Normal-Superconducting 
(SNS) junctions have been customarily interpreted using the theory j-of multiple 
Andreev reflections, following Klapwijk, Blonder, and Tinkham (KBT)EJ. 

KBT theory treats two conventional superconductors with a single s-wave BCS 
quasiparticle gap A, separated by a free electron metal. Electrons traversing the 
metal are Andreev reflected back as holes, gaining energy increments eV at each 
traversal (as depicted in Fig. [j]). Peaks in the differential resistance appear at 
voltages 2A/ne, and are due to the (E — A) -1 / 2 singularity in the quasiparti- 
cles' density of states. However, in cuprate SNS junctions, such as YBa2Cu3C>6.6 - 
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Fig. 1. KBT Theory. Differential resistance peaks of n = 6 (left diagram), and n = 5 (right 
diagram), involve a cascade of n Andreev reflected charges traversing the normal metal. Singular 
dissipation is due to emission of quasiparticles above the s-wave gap. Filled (empty) circles denote 
electrons (holes) in the normal barrier. 

YBa2Cu2.55Feo.45Oy -YBa 2 Cu306.6 examined by Nesher and KorerJij , application 
of KBT theory is problematic. A naive fit to KBT expression faces the two gaps 
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puzzle, i.e. an "AndreeAj-gap" is of order A « 16meV, while the tunneling gap is 
about three times larger!l3, and scales differently with T c . Without perfect aligp-i 
ment of the interfaces, it is hard to understand the observed sharpness of peakalj 
since the d-wave gap is modulated at different directions. Moreover, the barrier is 
by no means a "normal" metal devoid of interactions: it is an underdoped cuprate 
with antifcrromagnetic .correlations and strong pairing interactions as evidenced by 
a large proximity effectO. 

In the following sections we review an alternative explanation for the differential 
resistance peaks serieso, which takes into account the strong correlations in, the 
pseudogap regime. Our analysis resolves Duetscher's two energy scales puzzlca. 



5. SNS Junctions: pSO(5) Theory 

We consider a junction, where the barrier (N) has no superconducting or magnetic 



order (bl 



0. 



0. We derive on general grounds the form of the effective 



tunneling Hamiltonian between superconductors as follows. 

An integration of the barrier's charged bosons bh out of the path integral results 
in an effective action S tun which couples the charges of the two superconductors . 
S tun [bh L , bh R , b a ] explicitly depends on the hole pairs bosons on the left and right 
interfaces, and on the magnons in the barrier. By charge conservation, an expansion 

of S tun as a power series leaves only terms with equal number of bus and frt's. By 
spin conservation, the magnon terms are singlets and hence at least bilinear in n a . 
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Fig. 2. pSO(5) theory for Andreev peaks in cuprate SNS junctions. Three hole pairs 
co-tunneling from left to right, generate a pair of magnons. At the antiferromagnetic resonance 
threshold QeV = 2A S , this process contributes to the n = 3 peak of the differential resistance. 
The diagram contains lowest order contributions of hole pairs-magnon interactions to the tunneling 
vertex T3. 

This expansion leads to a series of tunneling terms. For the Andreev peaks we 
retain only the leading order terms (in b^ , b) which are 
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A — V^ T tf tf k h 

,/i ™ ~~ 2-^1 ±nU h L ,l- ■ - u h L ,n u h R n+l- ■ - u h R ,2n 

yi...j/2,!,x,x' 

x^n«(x)n Q (x')J (14) 

A^ describes a simultaneous tunneling of n hole pairs from the left to the right 
superconductor, coupled to a magnon pair excitation. T n is the tunneling vertex 
function, which depends on the bosons positions. 

The energy transfer mechanism is depicted diagrammatically in Fig.pl We do 
not compute T„'s which depend on the details of the barrier and the interfaces. A 
"good" N barrier is defined to have sizeable T n , if multiple pair tunneling terms 
are to be observed^-. This requires a thin barrier with slowly decaying spin and 
charge correlationsEH. It is important to note that multiple pair tunneling, i.e. the 
differential resistance peaks at n > 1, depends on strong anharmonic interactions 
between the hole pairs and magnons. These interactions are an essential part of 
the pSO(5) theory as modelled by H int in Eq.(g). 

The junction's conductance is calculated in the standard fashior£3: the bias 
voltage V transforms the left bosons bu L ~ * e t2eVt bh L , which yields time dependent 
operators A n (t). The current is calculated by second order perturbation theory in 

n tun-mag yielding 



AT'M = i dte^([A n (t),Ai}) (15) 



For singular contributions I sm s j W e ignore superconducting condensate fluctuations 
b\ — (bl), which have a smooth spectrum. Similarly, we ignore the frequency de- 
pendence of T n (w). Setting b^ -> (b{) and b[ -» e l2eVt (b[) leads to 

I«ng = \J 2 ne £ <&i) 4 "|T„[q]| 2 

n \qx\<Tr/d,\q y \<Tr/W 

xQj2S(q,iuj + 2neV + iO+)S(-q,iw) (16) 

UJ 

where the barrier dimensions are dxW (see Fig.gJ), and y] is a Matsubara sum. 
For a nearly antiferromagnetic "N" barrier, T„(x — x') in (114]) decays slowly with 
the distance between magnons. Thus for a narrow barrier a~<< W, the magnons 
are excited at q y rs 0, and the momentum sum reduces to a one dimensional sum 
over q x . At zero temperature we obtain 
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The last expression emphasizes the singular form of I sm s{y i A s ) at the peaks. For 
a large background conductance dl/dV » dl smg /dV, the inverse square root 
singularities in I sm s create peaks in the differential resistance dV/dl at voltages 



V n = A„/(ne), 



1,2, 



Q„ = 2ne 



(18) 



where Q n is the excess tunneling charge below the n-th peak. Note that Q n changes 
in increments of 2e. The differential resistance peak series is depicted in Fig. 0, for 
weak broadening of the singularities and an arbitrary set of coefficients t n . 
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Fig. 3. pSO(5) Andreev peaks. Eq. [n| is plotted for a choice of t n /A 3 J 2 = 2-™10~ 4 , n < 5, 
and a background conductance of unity. Below the n-th peak, the excess tunneling charge is 2ne, 
rather than BTK's ne. 



6. Discussion, and Proposed Experiment 

We have seen that magnon pair creation induces peaks in the differential resistance 
which arc similar in appearance to the Andreev peaks of the KBT mechanism. The 
crucial difference is that here the singular dissipative process does not involve Cooper 
pair breaking, but low energy antiferromagnetic excitations. In the KBT mecha- 
nism, a single sharp gap-like feature can be obtained in a d wave superconductor 
only by precise alignment of the a-b axes of the two superconductors. 

Here, one only requires the junction to be flat in the transverse direction, such 
that q y is conserved and the charge pairs are coupled mostly to the one dimensional 
singularity of the magnon density of states. This requirement is less stringent for 
weakly dispersive magnons near the resonance. 

In KBT theory for two identical superconductors, the peaks appear at voltages 
V^ BT = 2A/(ne), n = 1,2,... which are the upper threshold for tunneling of 
charges Q n = ne. Thus, KBT allows both even and odd number of electron charges 
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to participate in the multiple Andreev reflection process, as depicted in Fig. p], 
while the pS0(5) theory expects only pair charges Q n = 2ne. 

Observation of Andreev reflection enhanced shot noise S(V) has been reported 
by Dieleman et. a/J13 in a conventional SNS junction. They have measured the 
tunneling charge via the relatione!] S — 2Q n I(V n ). The increment of charge at the 
first Andreev peak at 2A was clearly seen to be of magnitude e. 

We propose that a similar measurement in YBCO junctions could provide a 
decisive discrimination between the processes of Fig.[j] and FigJ|. The goal is to 
measure the charge increments Q n — Q n -i at any peak position V n , n = 1,2,... 
and see whether they are of magnitude 2e rather than e. The measurement would 
probably involve a careful subtraction of the large but smooth background quasipar- 
ticle contribution to the current and the noise spectrum. We eagerly look forward 
to results of such experiments. 
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